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Abstract. Markovian regime decoherence effects in quantum computers 
are studied in terms of the fidelity for the situation where the number of qubits 
N becomes large. A general expression giving the decoherence time scale in 
terms of Markovian relaxation elements and expectation values of products of 
system fluctuation operators is obtained, which could also be applied to study 
decoherence in other macroscopic systems such as Bosc condensates and super- 
conductors. A standard circuit model quantum computer involving three-state 
lambda system ionic qubits is considered, with qubits localised around well- 
separated positions via trapping potentials. The centre of mass vibrations of 
the qubits act as a reservoir. Coherent one and two qubit gating processes 
are controlled by time dependent localised classical electromagnetic fields that 
address specific qubits, the two qubit gating processes being facilitated by a 
cavity mode ancilla, which permits state interchange between qubits. With a 
suitable choice of parameters, it is found that the decoherence time can be made 
essentially independent of N. 



1 Introduction 



The topic of quantum computation has developed enormously since the theoret- 
ical work of Feynman £Q and Deutsch [2] in the 1980s. Along with topics such 
as quantum teleportation, quantum cryptography and quantum measurement 
theory, quantum computation is one part of the expanding field of quantum 
information science jSj- Much of the current interest has been stimulated by po- 
tential applications of quantum computers in situations where they are expected 
to out-perform classical computers, such as in quantum algorithms for factor- 
ing large numbers 0] and searching large data bases [S], or in simulating the 
behavior of quantum systems 1 . The implementation of quantum computers 
has mainly focused on the standard quantum circuit model. Here quantum in- 
formation is stored as entangled states of an array of two-state systems (qubits) 
which are initially prepared with all qubits in one state, the algorithm is then 
implemented as a sequence of unitary operations involving one or two qubits at 
a time, and the result for the computation is provided by a measurement on the 
final state of the quantum computer. Other approaches have also been proposed, 
including adiabatic quantum computation 0, continuous variable, topological 



or holonomic quantum computation [7], [5], |§] and quantum computation based 
on a sequence of projective measurements instead of unitary quantum gates |10|. 
[TT] - including those using cluster states A variety of physical realisations 

of the quantum circuit model have been proposed, differing in the nature of the 
qubit system (e.g , hyperfine ground states and metastable excited states in ions 
and neutral atoms, one photon states in optical systems with two polarisation 
modes, ...), the way the qubits are prepared in the initial state (e. g., optical 
pumping - ionic and neutral atom qubits), the method used for the gating pro- 
cess (e. g., two photon- resonant laser fields - one qubit gates in ions and atoms, 
one photon resonant laser fields combined with photons in high Q optical cav- 
ities - two qubit gates for atoms in cavity QED based systems, cold collisions 
between qubits - two qubit gates for neutral atoms, ..) and the technique used 
to detect qubit states (e. g., quantum jump techniques for ionic qubits). Small 
scale demonstrations of quantum computation have already been achieved, such 
as implementing the Deutsch-Josza algorithm 13, and factoring small numbers 
|14|. Quantum computer architectures, which specify the components and how 
they are integrated within a coherent plan, provide useful guides for the long 
term research program needed to develop practical quantum computers, and al- 
low for developments of the basic architecture to overcome problems revealed as 
the research program is implemented. For example, ion trap quantum comput- 
ers have significantly evolved from the original Cirac-Zoller single trap proposal 
|15| to an architecture involving qubits shuttling between a memory region and 
a processor region in order to build a larger scale system |16|. The implemen- 
tation of quantum computers via quantum optical systems (atoms, ions and 
photons) is one of the more promising routes to follow , ^H] and the present 
paper deals with such a system. A set of criteria for the successful realisation 
of circuit model quantum computers has been formulated by DiVincenzo |T5] . 
and comprehensive surveys of current knowledge of the subject are contained in 
recent reviews and textbooks 

The idealised unitary evolution of the circuit model quantum computer does 
not occur in reality because the quantum computer interacts with the environ- 
ment, both during gating processes and when no gating is taking place. The 
loss of unitarity is referred to as decoherence, and can be quantified in terms 
of the fidelity. This specifies how close the actual behavior of the qubit system 
density operator is to its idealised behavior when only coherent gating processes 
take place and system-reservoir interactions are switched off. The time for the 
fidelity to change significantly from unity defines the decoherence time scale. 
Other measures of decoherence time scales have also been proposed [21] • Deco- 
herence is also of interest in other macroscopic systems such as Bose condensates 
and superconductors, and its general effect in macroscopic systems is important 
in quantum measurement theory and to understanding how classical behavior 
emerges |22| . |23| . In the case of quantum computers, one of the DiVincenzo cri- 
teria for effective quantum computation is that decoherence time scales must be 
much larger than gating time scales and the time intervals when gating processes 
are absent, and hence a determination of the decoherence time scale is important 
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in evaluating any specific quantum computer architecture. Decoherence is the 
enemy of quantum computation and a number of methods have been proposed 
for combating its effects. These include active (or error correcting) methods - 
quantum error correction [22], |2S], [2H|, dynamical decoupling techniques |27) . 
and passive (or error avoiding) methods - decoherence- free subspaces [2B]: EH]- 
PU) . [3*T] . topological or holonomic quantum computing [7J, [SI, 0- Combined 
methods, such as decoherence-free subspaces plus dynamical decoupling have 
also been proposed Quantum error correction implies a large overhead of 
redundant qubits and dynamical decoupling involves control pulses with time 
scales much shorter than reservoir correlation times, so in general terms it is 
desirable to implement error avoiding methods and use error correcting as a 
back up. It is well known that if the loss of fidelity in each gating process is kept 
below a certain threshold (estimated as being between 10 -3 and 10 -4 - see |33| . 
|2()| ) then fault-tolerant quantum computing is possible using error correcting 
codes [S3, [SSI- Error avoiding methods will be the focus of the present paper, 
with the aim of keeping the fidelity loss below such a threshold. The scalability 
of a quantum computer architecture is another important DiVincenzo criterion. 
The large qubit case is important for implementing quantum computers in sit- 
uations where they are expected to out-perform classical computers, such as 
factoring large numbers and searching large data bases 0. In this paper, the 
primary aim will be to study decoherence effects in quantum computers for the 
situation where the number of qubits N becomes large. 

In general |36| . the temporal behavior of a quantum system coupled to a 
zero temperature reservoir is of three distinct types, depending on the time 
regime: (a) Quadratic behavior at short times (t <C r c ) (b) Exponential decay 
at intermediate times (rd w f > r c ) and (c) Power law behavior at long times 
(t » Td). Here r c is the reservoir correlation time, Td the system decay time. At 
non-zero temperature T, a thermal time scale tj, = h/ksT is also involved. The 
short time regime is associated with the Quantum Zeno effect, the intermediate 
time regime involves Markovian decay and the long time regime behavior is due 
to a lower bound in the energy spectrum. 

For quantum computers, the physical relevance of the short time regime is 
not clear, since (apart from architectures making use of dynamical decoupling 
methods for error correction purposes |27| 'l most feasible measurements and 
gating processes are likely to require time scales much longer than the reser- 
voir correlation time. As Markovian theories for the intermediate time regime 
indicate that decoherence times decrease rapidly as the number of qubits N 
increases, the decoherence time may finally become comparable to the reservoir 
correlation time. Apart from creating consistency problems for the theory, this 
could place a limit on quantum computer size. The relationship between the 
short and intermediate time regimes has been studied for a simple case of N two- 
state systems, all initially in the lower state and coupled to the electromagnetic 
(EM) field in the vacuum state. A perturbation treatment correct to second 
order in the coupling constants [37] shows that the fidelity initially decreases 
quadratically for times less than t c , then reaches a minimum and eventually 
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returns to unity. The Markoff theory fidelity equals one at all times. The effect 
is due to the non RWA terms in the system-reservoir coupling, and may merely 
reflect the artificial nature of the uncorrelated system-reservoir initial state. 
Alternative treatments of the short time regime based on methods preserving 
unitarity ((SSj, [SH]) may overcome these difficulties for studying systems with 
very large N. A preliminary study |4U) on a standard quantum computer model 
with two-state qubits obtained a quadratic behavior of the fidelity in the short 
time non-Markovian regime. 

The present work deals with the intermediate time regime. The interme- 
diate time regime is more physically relevant for reasons discussed above (the 
quantum computer model we study does not involve the use of dynamical de- 
coupling), and here Markoff theory can be used. For internal consistency, the 
decoherence time scale td must be long compared to r c . Markovian expressions 
for the intermediate time behavior of the rate of change of fidelity are obtained 
for the general case where the qubit system (including any ancilla) are in a pure 
state and the reservoir (which may involve several components) is in a thermal 
state. The initial rate of change of the fidelity defines the decoherence rate and 
its inverse is the decoherence time scale td. This result gives the decoherence 
time scale in terms of Markovian relaxation elements and expectation values 
of products of fluctuation operators for the decohering quantum system. The 
expression is quite general and may have applications for decoherence in other 
macroscopic systems, such as Bose condensates or superconductors. For the 
quantum computer model studied, the characteristic decoherence time scale is 
evaluated at finite temperature for specific qubit states (such as Hadamard and 
GHZ states) in the situation of no gating processes occuring (memory decoher- 
ence). Decoherence effects on quantum computers due to one and two qubit 
gating processes are studied for zero temperature (gating decoherence). The 
zero temperature case should be most favourable for long decoherence times, 
and short decoherence times even at zero temperature during gating processes 
would be ominous for implementing large scale quantum computers. 

As mentioned previously, architectures involving ionic qubits are amongst 
the most promising for possible implementations of circuit model quantum com- 
puters \T7\, [201 • However, directly scaling up the original Cirac-Zoller 
model ,15j to large numbers of qubits is difficult not only because it is hard to 
trap large numbers of ions in a single trap, but also because using a collective 
vibrational mode as an ancilla for two qubit gates becomes impossible when the 
mode density becomes very large. An extensive study of decoherence effects as 
qubit numbers are increased in Cirac-Zoller type quantum computers where a 
vibrational mode is used for gating purposes, has been carried out by Plenio and 
Knight 021 f° r both two state and three state lambda system qubits, and 
applied to real ions. Limitations on the size of quantum computers (and there- 
fore on the size of numbers factorisable via the Shor algorithm 0]) based on this 
model were found. Even allowing for error correction this limit was quite small. 
A solution to the problem of scaling up qubit numbers is being developed using 
shuttled qubits [16| . and a different approach involving a 2D array of trapped 
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ions with a control qubit moved above its target target qubit to carry out a two 
qubit gate via a collisional process has also been suggested 03] > EH- The 2D 
arrays might be based on elliptical rather than linear ion traps 05J- Another 
approach not involving ion shuttling has also been proposed, with two qubit 
gating via Raman laser pulses applied to the pair of qubits [4"B], Architectures 
involving high Q cavities are also attractive ^7], HHIi an d it is worthwhile 
to try to combine these with ionic qubits. Cavity QED systems involve using 
an optical (or microwave) cavity mode photon as an ancilla to enable two qubit 
gating to occur for atomic or ionic qubits, and were in fact amongst the earliest 
proposals for quantum computers 07|, |3S], 0H|. High-fincssc optical cavities 
have now been developed [3U| which allow a large number of atomic qubits to be 
enclosed under strong coupling conditions, so that coupling of the cavity mode 
to the qubit system is faster than both cavity photon loss and qubit decay 
via spontaneous emission. In such circumstances, two qubit gating via cavity 
photons can occur much faster than decoherence due to the loss processes. How- 
ever, it is difficult to create strong coupling conditions for ionic qubits, since the 
small optical cavity volume can interfere with the ion trap and the presence of 
electrodes can interfere with the cavity mode though single ion trapping 
inside a cavity under weak coupling regime conditions has been realised |52| . 

. There have been several proposals for scalable quantum computers involv- 
ing ionic qubits in optical cavities under weak coupling conditions |3I], [S^, 
|56|. |51|. involving probabilistic entanglement protocols or using dissipation to 
confine the evolution in decoherence-free subspaces. The optical cavity may be 
arranged with its axis perpendicular to the array of ions, so that two qubits 
at a time are in the cavity. Strong coupling regime proposals have also been 
formulated [52! ■ In the present paper, we consider a combined ionic qubit and 
high Q cavity quantum computer architecture with a large number of qubits 
in the cavity. Our model is somewhat similar to the two qubit case considered 
by Tregenna, Beige and Knight |56| . but now we consider a large number of 
qubits and also allow for their vibrational motion rather than treat them as 
stationary. As the emphasis of the present work is to examine the effects of 
decoherence on the system considered by Tregenna et al |5j)] as qubit numbers 
are increased, the same parameters as in their work will be used, rather than 
those for real ions. The issue of developing a theory for real ions is discussed in 
the last section of the paper. Memory decoherence for the Cirac-Zoller model 
[T5] due to the effects of vibrational motion has also been studied by Garg 
Note that in the present case (unlike in the work of Plenio and Knight |42|). 
the vibrational modes act as a reservoir rather than as an ancilla to facilitate 
gating processes. 

A standard model involving N ionic qubits is considered, the overall ar- 
chitecture being illustrated in figure 1. Each qubit is in a three-state lambda 
system rather than a two-state system as previously treated in the short 
time regime |4(J|. The qubit states are the two lower states 0, 1. The quantum 
computer system also includes a high Q cavity mode, which is coupled to the 
qubits and acts as an ancilla. Lambda systems, as well as facilitating Raman 



5 



gating processes, should result in qubits that are less vulnerable to spontaneous 
emission (SE) based decoherence, the upper state 2 only being occupied during 
gating processes. Reducing SE by having qubits located in a high Q cavity is also 
desirable, and the cavity mode also facilitates two qubit processes. Our model 
involves N ionic qubits localised around well-separated positions via trapping 
potentials, and the centre of mass vibrational motions of the qubits are now 
treated. Coherent one and two qubit gating processes are controlled by time 
dependent localised classical EM fields that address specific qubits. The one 
qubit gating process involves weak two-photon resonant Raman gating fields, 
well detuned from one photon resonance. The two qubit gating processes are 
facilitated by the cavity mode ancilla, which permits state interchange between 
qubits. For the two qubit gating process, resonant gating fields coupled to the 
2-1 transition for the control qubit, and coupled to the 2-0 transition for the 
target qubit. The cavity mode is resonant with the 2-1 transition, but uncou- 
pled to the well-detuned 2-0 transition (as in [EE]). Two qubit gating processes 
take place in decoherence- free subspaces [EI], [EE]- In our model, the reservoir 
(or bath) has three constituents. The three-state qubits are coupled to a bath 
of spontaneous emission modes, and the cavity mode is coupled to a bath of 
cavity decay modes. For large N the numerous vibrational modes of the ionic 
qubits also act as a reservoir, with Lamb-Dicke coupling to the internal qubit 
coordinates, cavity mode and gating fields. Non-RWA couplings are included. 
All qubit interaction terms (electric dipole, Rontgen, diamagnetic, ionic cur- 
rent) with these three baths are examined and the important contributions to 
the decoherence rate found. The qubit-bath coupling is amplitude coupling via 
cr± optical coherence operators. In addition to these fundamental causes of de- 
coherence, technical shortcomings in the implementation of the computer model 
can also cause decoherence. For example, the trapping potential producing the 
array of qubit trap sites could be subject to fluctuations. Such trapping poten- 
tials could be provided by off resonant near-classical optical fields or by magnetic 
fields, and these could fluctuate. Decoherence effects due to fluctuations in these 
fields could have significant effects [EE] and should be evaluated . However, we 
will concentrate in this paper on fundamental causes of decoherence, especially 
the effects of qubit vibrations, and technical causes will be left to a later time. 

The plan of this paper is as follows. In section 2 we set out the Hamiltonian 
for the quantum computer model and derive our general expression for the 
decoherence time scale by applying Markovian evolution theory. In section 3 
the decoherence time scales are evaluated for the no-gating, one qubit gating 
and two qubit gating cases. A summary of the main results of the paper is 
presented in section 4 along with a discussion of extensions of the theory for 
real ions. 



2 Theory 
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2.1 Hamiltonian 



The total Hamiltonian for the system can be written as 

H = H s + H c + H B + V s + Vj (1) 
where the Hamiltonian for the qubit system and cavity mode ancilla is 

H s = J2^i a + tuJ b b% (2) 

i a 

the Hamiltonian for the collective vibrational motions of the qubit system is 

ia ijot/3 

= hv K A ] K A K (4) 

K 

and the Hamiltonian for the bath of spontaneous emission and cavity decay 
modes is 

H B = Y^ hw k a\a k + Y H k b\b k . (5) 

fc k 

The coherent coupling Hamiltonian for gating processes in the qubit system 

is 

V S = WUa+nlaX*++*-) 
i;a=0,l 

+ E Hg t ab + g* a tf)(^ + a™) (6) 

j;a=0,l 

The qubits are coupled to both classical fields and the cavity mode. Each qubit 
is addressed by localised classical EM fields to facilitate 1 qubit and 2 qubit 
gating. For 1 qubit gating Q ia for ith (gated) qubit are two photon resonant 
Raman fields strongly detuned from 0-2 and 1-2 transitions. For 2 qubit gating 
flu is resonant with the 1-2 transition for the ith (control) qubit and f^o is 
resonant with the 0-2 transition for the jth (target) qubit The cavity 

frequency ujf, is resonant with the 1-2 optical transition [SHj. In the model of 
Tregenna et al the cavity mode is only coupled to the 1-2 optical transition 
9io = 

Finally, the interaction of the qubit system and ancilla with the bath and 
centre of mass vibrations is given by 

Vr = E ^+"-)E^°* + sM) 

i;a— 0,1 k 

+htf[J2(vkb k + w k b\) + J2 tx(A - A K )] + HC 

k iK 
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+ ]T ft(<7* + ° + a™)bJ2PK(A K + rf K ) + HC 

i;a=0,l K 

+ %T + ^)E( 8 UeS(4 + 4) (7) 

i;a=0,l if 

Terms include electric dipole coupling of qubits to SE modes, quasi-mode cou- 
pling of cavity mode to decay modes, Lamb-Dickc coupling of qubits, cavity 
mode, gating field to CM modes. Rontgen and diamagnetic terms are not in- 
cluded as their effects were shown to be small. 

The centre of mass (CM) displacement of the qubits is related to the collec- 
tive vibrational modes of the qubits via 

Sr ia = T, S ^ k \[^^( Ak + A k) ( 8 ) 
£^f%3;K - mv 2 K S ia , K (9) 

30 

where the unitary real matrix S relates qubit CM displacements Sr ia (a — 
x, y, z) to vibrational normal coordinates. 

The qubit, cavity mode, centre of mass motion, bath modes operators are 



a + 


= (|2)(a|) 4 


a™ = (|a)<2|) i 


a ab 


= (\a)(b\)i 


a 7^ b a = 0, 1 


a% aa 


= (\a)(a\)i 


a = 0,1,2 


[6,6t] 


= 1 [ak,a\ 


] = s H 




= S ki [Ak: 


4] = Skl 



a = 0, 1 

(10) 
(11) 

These include qubit optical, Zeeman (or hyperfine) coherences and population 
operators, as well as bosonic annihilation, creation operators for the cavity, 
spontaneous emission, cavity decay and the CM vibrational modes. 

The Hamiltonians involve certain coupling constants defined as follows: 

= -»£W 2^fiy ^ 2a ' Uc ~) ac ex Pi ik c ■ no - u c t) (12) 

r. * 



3ia = -i\j^^-(d 2 a-u b )exp(ik b -r l0 ) (13) 



9k = 2 ^ k h y ( d 2a ■ u k ) exp(ik ■ r l0 ) (14) 



/ LU b 


1 H 


2e hV b \ 


' 2mvK 



Pk = \ 7, — rrr\/o ( d 2a • u b )(k b ■ S iK ) x 
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x exp(ik b ■ r i0 ) (15) 



xa c expi(fc c • r i0 - oj c t) (16) 

4 = ieT 7liw/S (fcb -^ )cxp(ifcb - rjo) (17) 

where a = 0, 1 unless stated otherwise. The cavity mode volume is V&, the SE 
mode volume V. Each ion has charge e-r- 



2.2 Dynamics and Decoherence Time 

The total density operator for the qubits, ancilla and reservoirs satisfies the 
Liouville-von Neumann equation 

dW 

l h— = [H,W] (18) 
H = Hs + Vs + Hr + Vj (19) 



The initial condition is given by 

W(0) - P s (0)p R (0) (20) 

and represents an uncorrelated state for qubits and reservoirs, the qubits initially 
being in a pure state \ip s ) and p s (0) — \il> s )(ip s \, whilst the reservoirs are in 
thermal states. 

The reduced density operator for qubits and ancilla is defined as 

Ps - Tr R W (21) 

and its general evolution allows for both coherent coupling and reservoir inter- 
actions. 

For coherent evolution due to Vs only, the reduced density operator for 
qubits and ancilla would satisfy 

lh ^ = [Hs + V s ,p S0 ] (22) 

where the same initial condition Pg (0) = pg(0) can be chosen as for the general 
evolution. 

The fidelity is defined by 

F = Trsp so p s (23) 
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and compares the actual and idealised evolution of the qubit system. 

The decoherence timescale is defined as the time scale over which actual 
quantum computer evolution (non-unitary) differs significantly from idealised 
coherent evolution (unitary) . The decoherence time scale will be defined via the 
time dependence of the fidelity. It is related to certain basic time scales due 
to qubit-environment coupling. These are: (a) the reservoir correlation time r c 
- for EM field SE modes t c ~ 10~ 17 s (b) Markovian relaxation times - T\ for 
populations, Ti for coherences - for EM field SE modes T\,i ~ 10~ 8 s (c) the 
thermalisation time r b — h/ksT - At lpK r b ~ 10~ 5 s. The decoherence time 
scale will also depend on factors such as qubit system state, the numbers of 
qubits and the reservoirs involved. 

Markovian evolution occurs for t » r c and the reduced density operator 
satisfies a master equation 

[H s + V s ,p s ]+Lp s (24) 
~iJ2A ab [S a Sl Ps } 

ab 

+ E^a b {[Sl Ps S a } + [Stp s ,S a }}. (25) 

ab 

The master equation includes the Liouvillc supcroperator term involving Marko- 
vian relaxation T a b and frequency shift A a b matrices and system operators S a , 
where the system-reservoir interaction Vi is sum of products of system and 
reservoir operators 

Vj = J2SaRa (26) 
a 

[H s ,S a ] = huj a S a . (27) 

The Markovian evolution requires the reservoir correlation functions (R a (t)R b (t— 
r)t) approach zero for r » r c , and the Markovian relaxation T a b = and 
frequency shift A a b = AJ a matrices are given via reservoir correlation functions 

as: 

C a b = J dT{R a {t)R b {t-T)^)e^>{-{iuj b + e)T) (28) 
o 

r a6 - (C o6 + C b * a )/2, A ab = (C ab - C* ba )/2i (29) 



In the Markovian intermediate time regime (t^ sj ( » t c ) the rate (dF/dt)o specifies 
the decoherence rate for the qubit system, and its inverse defines the decoher- 
ence time r d The decoherence time can be expressed in terms of the Markovian 
relaxation matrix elements and qubit system averages of products of fluctuation 



ih 



d Ps 
dt 

Lp s 
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operators 



t d = { at 



= 2(j2T ab (AS a ASl)s) (30) 

t»T c ^0 \ab J t> rc ^0 



Here the qubit and ancilla are in pure state \xs)> which may differ from ini- 
tial state \rp s ) due to gating processes. For t » t c ^ we take p so — p s — 
\Xs)(Xs\- In tne above AS a = S a - (S a ) s and (D) s = Tr s (np s ). The expres- 
sion for the decoherence time scale is quite general and may have applications 
for decoherence in other macroscopic systems, such as Bose condensates or su- 
perconductors or in the theory of quantum measurement. A special case of this 
result is given in Ref. |§Jj], which deals with Bose condensates. 

Expressions can be obtained for the loss of fidelity. During a time T much 
smaller than r jj , the change in fidelity if no gating is occuring is 

AF — ( ) no — gating -T (31) 
T D 

and if gating is taking place during a time AT the change in fidelity is 

AF = -{ — )g atmg .AT (32) 
td 

For idealised quantum computation we require A_F <C 1. 

In the short-time regime (t <ti t c ) the time dependent fidelity may be written 
in a power series |4l)| 



F(t) = !-(!)_(* ) + .. (33) 
Tl 2t 2 



Tl 

r, 2 



(34) 



5-5 - Tr BC (il>s\ViWW s )p B (0)p c (0) 

ZT 2 

-rr BC (^ 5 |y / (0)|^ 5 ) 2 p S (0)p c (0) (35) 
= Tr BC ((AV I (0) 2 ) s )p B (0)p c (0) (36) 

The qubit and ancilla are initially in a pure state \tpg) and Psi®) — \'4 ) s)(' i l ) s\- 

Here AVr(O) = Vr(0)-(Vr(Q)) s where (Q) s = Tr s {np s ). The times t 1 ,t 2 , ..specify 
characteristic decoherence times for the qubit and ancilla system, their inverses 
defining decoherence rates. For the uncorrelated initial state W(0), only T2 is 
involved in specifying the short time decoherence. The decoherence time scale 
T2 depends on qubit and reservoir averages of the fluctuation in the system- 
reservoir interaction operator squared. The quadratic time dependence of the 
fidelity is characteristic of the quantum Zeno effect. However, results obtained 
for T2 are due to non RWA terms in the system-reservoir coupling, and may 
merely reflect the artificial nature of the uncorrelated system-reservoir initial 
state. 
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3 Results 



3.1 Case of No Gating 

For the situation where no gating processes are occuring, there is no upper state 
1 2) amplitude, the cavity mode is in a no photon state \0)a and the qubit-ancilla 
state \xs) is given by I^q)^)^, where qubit system state is |0q). This situation 
corresponds to states produced after idealised coherent gating processes. To ex- 
amine unfavourable scenarios, the reservoir temperature T is assumed non-zero 
and the qubits and the cavity is assumed to have a low Q, so that spontaneous 
emission decay leads to a larger decoherence rate than would otherwise be the 
case. Spontaneous emission is the dominent decoherence process, and only its 
contribution is shown. 

In the Markovian intermediate time regime {jd ~ t > r c ) the decoherence 
time is given by 

— =exp(-|^)EVr^cos0 Qb E« 6 ) (37) 

T D KB J- ab i 

The expression involves Zeeman (or hyperfine) coherences for ith qubit (cr ab ) (a ^ 
b) (a, b = 0, 1) and populations for zth qubit (<r aa ) (a — 0, 1). The optical tran- 
sition frequencies are t02 a ~ u>o (a — 0, 1), the spontaneous emission decay rates 
for 2-a transition are T a and the angle between dipole matrix elements d,2 a , d-ib 

is e ab . 

For the case of the Hadamard state (uncorrelated), the qubit state is \4>q) = 
I^q)») where i th qubit state vector is \<t>Q)i = (|0) ? ; + |l)i)/2 1 / 2 , we find that 

i 

± = l N eM-^)EVraT b co S ab (38) 
td 2 k B L ab 

and note that theecoherence time for Hadamard state can be infinite for the 
lambda qubit system if choose e^o + <^2i = 0. 

On the other hand, for the case of the GHZ state (correlated), where the 
qubit system state vector is \4>q) = (|00..0) + |ll..l))/2 1 / 2 , we obtain 

- = ^«p(-^)Er. 09) 

TD 2 kB± a 

Here the decoherence time scale is still very long, due to the upper state Boltz- 
mann factor. With N «10 4 qubits and T wl0 8 s _1 , lo wlO^s" 1 , T w300 K, 
we find td ~10 19 s. 

Overall, spontaneous emission due to electric dipole coupling is dominent 
cause of decoherence. Other terms such as Lamb-Dicke coupling can be ignored, 
so the qubits can be treated as stationary. The decoherence time scales as 1/N 
for this intermediate time regime but it is very long. The present results for this 
case of memory decoherence are consistent with those of Garg 58 j; who also 
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found very low decoherence rates allowing for vibrational motion in Cirac-Zollcr 
type quantum computers |15| with large numbers of qubits. 

3.2 Case of One Qubit Gating 

For the one qubit gating process with two-photon resonant Raman gating fields 
detuned from one photon resonance, the upper state \2) amplitude for ith 
(gated) qubit becomes slightly non-zero, though the cavity mode remains in 
the zero photon state |0)a \Xs) gi yen by \4>q)\0)a, with the qubit state \4>q) 
(a,j — 0, 1) given as 

{a} {a,} 

Here {a} = {ai, .., a/v}> { a i} = {ai, ■, it-i, Oj+i, ■■}• However spontaneous 
emission decay is small in the high Q cavity. The reservoir temperature is 
assumed zero. 

In the Markovian intermediate time regime {jd ~ t 3> r c ) the decoherence 
time is given by 

= 2{ J2( a ab)^ka-C;kb-C 

T D k^i ab 

+E(Kb)-(v M )(°tma-,.b- 

ab 

+E((4 2 )^-K a )R 6 ))r ia+;ft+ 

ab 

+ Y,( lT ab) T *ia-C;ib-C + Yl (vh^ia+C; ia+C } (41) 
ab a 

The decoherence time involves Markovian relaxation elements and state depen- 
dent quantities for the qubit system. For the gated qubit terms, optical co- 
herences (ct±), Zeeman (or hyperfine) coherences (cr ab ), upper state population 
((T22 ) an d lower state populations (a aa ) are involed. For the non-gated qubits 
[k 7^ i), Zeeman (or hyperfine) coherences (c^ & ) and lower state populations 
(°"oo) are present. 

Expressions for the zero temperature relaxation matrix elements are as fol- 
lows: 

T ka -c-,kb-c = lv 2 9k a gib'^T (k^i,k = i) (42) 

^ia-,ib- = l -V 2 (^* ib -n ia n ib )/A° ci (43) 
Fia+;ib+ = 7;V^a^bCOs8 a b 
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+~ r] 2 (rL* a n* b -n ia n ib )/A ci (44) 

T ia+C]ib+ C = ir?g la g* b {l ~ 5ab){ ~ l) (45) 

^max 

In these formulae the Lamb-Dicke (LD) parameter is 77, the gating EM field 
Rabi frequencies are f2i a (a = 0, 1), the one photon detunings for the Raman 
gating field are A° ci and the one photon Rabi frequencies for the cavity mode are 
9ka (a = 0, 1). The Zeeman (or hyperfme) transition frequencies are uj a b (a ^ 
b) (a,b — 0,1). Vibrational frequencies range from zero up to ^ ma x, and the 
vibrational modes have zero phonons at absolute zero. In the case studied the 
cavity mode is resonant with the 2-1 transition = LO21 ~ u->o, but both one 
photon Rabi frequencies gka are assumed non-zero. Approximations based on 
A°j > f max , T and wi > ^ ma x have been used. 

The populations and coherences are found for large one photon detuning, 
assuming O^o = ^iiexpiA^ = £l(t), where the common amplitude f2 has a 
maximum fi m and a width AT. For a gated qubit initially in state |0) (see 
Vitanov et al [01]) we find that 





= Ko>* = 


i sin 


°oo) 


= cos 2 9, 


Ki; 


«°> 


= (<T M )* = 


n 

_ A' 




= = 


n 

~A' 


P22) 


= n 2 /A 2 





9. cos 9. expiA^ (46) 
= sin 2 9 (47) 

cxp(-i6i).cos6» (48) 

exp(— i9).ism9.exp iAtfi (49) 

(50) 

In these expressions the quantity 6(t) is defined by the integral of the two photon 
Rabi frequency SIr — f2 2 /A 

9{t) = J dt'n(t') 2 /A (51) 

—00 

The gating time AT is given by 9 — tt/2 for t ~ AT, corresponding to the time 
the qubit takes to evolve into state |1) 



The parameters used in both the one and two qubit gating processes are set 
out in Table 1 



Table 1. Parameters used in gating processes 
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s.io 1 ^- 1 










^max 





|O ia | 


S.lO^s" 1 




3.10V 1 


9ka 


3.10V 1 


r a 


3.10V 1 


V 


6.1Q- 2 



V 

x cav 


3.10V 1 


Q 


1.10'' 



These correspond to optical and hypcrfinc transitions with a high Q optical 
cavity in the medium coupling regime. As discussed previously, these param- 
eters are chosen to be the same as in the work of Tregenna et al so that 
their model can be compared to the present one where decoherence effects due 
to vibrational motion is allowed for as qubit numbers increase. The cavity cou- 
pling constant and cavity decay rate are made equal, and large compared to 
spontaneous emission (SE) rate, and the gating field one photon Rabi frequency 
is large compared to the SE decay rate Two cases are studied, corresponding to 
the gating field one photon Rabi frequency being (i) small compared to (ii) the 
same as the cavity coupling constant and cavity decay rate. Case (i) and with 
zero gko applies in Tregenna et al |56| . The maximum vibration frequency and 
Lamb-Dicke parameter are calculated for a 3D Ca + lattice with lattice spacing 
3fi. A standard approach to the theory of lattice vibrations is used [HI], in 
which the vibrational potential energy is given by a quadratic form of the small 
displacements of the qubits from the lattice sites, as in Eq[3] The interaction 
between each pair of ionic qubits is obtained from electrostatics, from which 
the V^ 13 are obtained. For the cubic lattice case, the vibration frequencies are 
obtained from EqEl and are approximately proportional to the magnitude of 
the wave vector for each vibrational mode. Expressions for the quantities Si a -K 
are obtained from Eq|5| 

Using these parameters, the zero temperature relaxation elements can be 
calculated and the results are given in Table 2. 



Table 2. Relaxation elements in one qubit gating process 



Tfca-C; kb-C 


Tfa — ; ib— 


r^a+; ib-\- 


Tm+C; ib+C 


3.10-V 1 


LIOVV.IOV 1 


S.lOV 1 


4.10V 1 






l.lQ u s - \ 1.10V" 1 





The overall contributions to the decoherence rate can be obtaining by com- 
bining the qubit factor (for 9 ~ 7r/4, midway through process) and the relaxation 
matrix elements for the various terms (1 — 5) involved in the expression for the 
decoherence rate 1/td. The results are presented in Table 3. Here NG, G re- 
fer to non-gated and gated qubits respectively. The two gating field cases are: 
(i) = 3.10 6 s _1 (ii) £1 = 3.10 8 s _1 .The corresponding gating times are 10 _2 s, 
10~ 6 s (see Ea l52f) . which are long compared to the correlation time for the vi- 
brational modes reservoir of 10~ 7 s (given by 27r/i/ max ), and much longer than 
the correlations times associated with the SE or cavity decay modes reservoirs. 



Table 3. Contributions to decoherence rate in one qubit gating process 
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Term 


Relaxation 


Relaxation 


Qubit 


Overall 1/td 




Effect Due 


Rate 


Factor 


Contribution 


1, NG 


LD — Cavity 


S.lO^s- 1 


0.5 


iV.2.10- y s- i 


2,G 


LD — Gating 


1.10°s _1 


0.5 


(i) h.lQ- l s- 1 






i.io^- 1 


0.5 


(ii) 5.lO +i s- L 


3,G 


SE 


a.io^- 1 


lO" 8 


(*) a.io-^- 1 






a.io^- 1 


io- 4 


(ii) S.lO^s- 1 


4,G 


LD — Cavity 


S.lO^s" 1 


0.5 


i.io-V 1 


5,G 


LD — Cavity 


4.10 t> s- i 


IO" 8 


(i) 4.10- 2 s- i 






4.10 t> s- i 


io- 4 


(ii) 4.10+^s- 1 



The scaling of the decoherence rate with qubit number for tne case of one 
qubit Raman gating gives an overall decoherence time which is essentially in- 
dependent of N. This is because the terms for gated qubit i do not scale 
with N, whilst the non-gated qubit contributions (which scale with N) are 
negligible in comparison, even for N ~ 10 4 qubits. These features are clear 
from the results in Table 3. For both the small and larger gating held cases, 
the largest overall contribution to the decoherence rate arises from the terms 
associated with Lamb-Dicke coupling of the gated qubit with the gating field 
12(( a ab) ~ ( <T - ) ( a +))^ia-; »6— • The next largest overall contribution arises from 

ab 

the terms associated with Lamb-Dicke coupling of the gated qubit with the cav- 
ity mode J2( a 22)^ia+C; ia+c ■ The larger relaxation rate T ia+c -ia+c for the 

a 

latter is balanced by the smaller qubit factor (a l 22 ) associated with the upper 
state population. 

The overall fidelity loss can then be obtained. Since terms associated with 
Lamb-Dicke coupling with the gating field arc dominent, these terms combined 
with the gating time gives for the overall fidelity loss 

AF ~ 



We note that the fidelity loss only depends on Lamb-Dicke parameter. For 
rj = 6.10 -2 , the fidelity loss is 2.10 -3 . This is reasonably small though still 
somewhat large for fault-tolerant quantum computation - see . |2()| . 

3.3 Case of Two Qubit Gating 

For the two qubit gating process, with resonant gating fields coupled to the 2-1 
transition for ith (control) qubit, and coupled to the 2-0 transition for the jth 
(target) qubit, and the cavity mode resonant with 2-1 transition, but uncoupled 
to the 2-0 transition [55], the upper state |2) amplitude for the ith (control) and 
jth (target) gated qubits could be non-zero. The cavity mode is in zero photon 



1 rfOL 
2' A 

n 2 



7T A 

2 W 



(53) 
(54) 
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state |0)yi when one or both gated qubits i, j are in the upper state, but could 
be in one photon state \1)a when both are in lower states. Here \xs) given by 
I0q>|O)a + \4>q)\1)a, with qubit states \<jP£) 

I/q> = E^(K}; 2 A)IK}; 2 i 2 ^)+E c o(K};2jlK};2,) 

{aij} {a;} 

+ E ^({a,-}; 2 J .)|{a J .}; 2 .}+ £ C ({a})|{a}) (55) 

{aj} {a} 

I^q> = E C i(W)IW)- (56) 

{a} 

Here {a} = {ai, a 2 , .., a N }, {aj = {ai, ., a;_i, a i+ i, ..}, {a,} = {ai, ., aj-i, a j+1 , ..}, 
{a^} = {ai, ., aj_i, dj+i, .., fflj-i, flj+i, ••}, and (a,,aj = 0,1). However sponta- 
neous emission decay is small in the high Q cavity. The reservoir temperature 
is assumed zero. 

In the Markovian intermediate time regime (r^ - t » t c ) the decoherence 
time is given as the sum of non-gated and gated qubit contributions as 



(57) 



-H- +(-, 

tdJ \ t dJ ng \tdJg 
The non-gated (NG) qubits contribution involves 7 terms 

~ D J NG kjti,j ab 

k^i,j ab 

+ E E((<Tabb)r k a-,kb-C++Cc) 
k^i,j ab 

+ E T,{<7ab bbt ) r ka-C;kb-C 
k^i.j ab 

+ E E(^ t &>r fc a-C + ; fe6 -C+}. (58) 
k^i,j ab 

The gated qubits contribution consist of 120 terms. Some of the terms are 

-, \(23) 

— = 2E(K fc )-(0« b »r M - ;lb - (59) 

T D/ G ab 
1 \ (24) 

— = 2E((4; 2 J-WK fc ))r»-;i t - (60) 

T D/ G ab 
1 \ (34) 

— = 2E(«2;26>-<0^+»ria-;i6- (61) 
T D/ G ab 
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(35) 

= 2E(K h )-(^ a )K"»r,a-; Jh - (62) 

~D/q ab 
1 \(120) 

— = 2«6tfc> - <6+)<6))r c+ ;c+ (63) 

t dJ g 

Some of these are equivalent to those for one qubit gating, others involve dif- 
ferent expressions, since the gating fields are now resonant rather than having 
a large detuning. In addition, there are many new terms only present for two 
qubit gating. Again, the decoherence time involves Markovian relaxation ele- 
ments and state dependent quantities for qubit system. The terms for non-gated 
qubits (k ^ i,j) - involve Zeeman (or hyperfine) coherences (<J k b ) and lower 
state populations ((T k a ), as well as quantities also involving the cavity mode op- 
erators (a k ab b), (a k b tf), (a k ab btf), (a k b tfb), (a k aa b), <cr* &t), < ff * 6&t) and (v k aa tfb). 
For the terms involving gated qubits, optical coherences (cj. a ), Zeeman (or hy- 
perfine) coherences {<J g ab ), upper state populations (a^), lower state popula- 
tions {a 9 aa ) (g refers to i, j), as well as two qubit state transitions (cr?>2-ab)> 
(°a&-22)' (°"a2-2b) an< ^ (^fa-fta) are a ^ involved. The two qubit transitions are 
a ab-cd = (l a )( c l)s(l^)(^l)hi where g ^ h refers to i ^ j). In addition, there 
are quantities involving cavity mode operators also, similar to those for the 
non-gated qubits. 

For the non-gated qubits, all relaxation elements r fea _ ; fcb _ , T ka _ ; k b-c, r fco _ ; kb-C+, 
Tfca-C; kb-c and Tka+c-, kb+c are zero for the specific gating process in- 
volved. In particular, the last two are zero because the cavity mode is resonantly 
coupled to the 2-1 transition and uncoupled to the 2-0 transition, as may be 
seen from the following expressions: 

tka-C;kb-C = g 7 ? 9ka9kb— (64) 

t \ ■ 2 (! - S ab)(- 1 Y l ( n r\ 

Tka+C;kb+C = l"f] 9kag kb (65) 



For the gated qubits (i control, j target qubit) some of the relaxation ele- 
ments are 

^ia-;ib- = i SalSbl ^(Qnflii - fi*if2*i)/^max (66) 

Tia- jb- = i SaiSbo i] 2 (k cl ■ k c j)(SlnQ,jo - H^Cljo) / v max 
1 Xi j , sin a; 



x 



_ / dx — (67) 



Xij Q x 



Fj a —,ib- — (Fia—,jb-)* (68) 
I*.;.. = i SaoSbO T] 2 (fljoQjO ~ ^ * £1 * ) / f max (69) 

Tc+;C+ = -j^cav (70) 
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The quantity Xij — y3a / |r*io — fjo\ relates the qubit separation to the lat- 
tice size As stated above, the cavity mode is resonant with 2-1 transition 
(ujb = u>2i ~ wo), the cavity mode is uncoupled to the 2-0 transition (gko = 0), 
the control gating field is coupled to the 2-1 transition, the target gating field 
is coupled to the 2-0 transition and both gating fields are on resonance. Ap- 
proximations based on loq ^> loiq 3> v max are used in the derivations. The 
vibrational modes have zero phonons at absolute zero. 

Expressions for qubit populations, coherences and two qubit transitions 
could be obtained based on the work of Tregenna et al 1^01 ■ Gating is based on 
the use of decoherence-free subspaces. States for the gated qubits (|0j0j), |0,lj), 
|lj0j), \A) = ( | li2j) — |2^1j))/\/2) 5 with the cavity mode in zero photon 

state \0)a are not directly coupled to one photon cavity states (see (SB]) via the 
qubit-cavity interaction. A CNOT gate can thus be performed with negligible 
cavity mode excitation, thereby avoiding decoherence due to cavity mode decay. 
Their treatment assumes = \£ljo\ = O(t), where f2 has a maximum VL m 

and a width AT. 

The gating time is given by 

AT ~ (71) 

The parameters used for the two qubit gating case are the same as for one 
qubit gating, except in accordance with Tregenna et al, the gating field is weak, 
and only Sl m = 3.10 6 s _1 is used (see Table 1). The corresponding gating time is 
2.10 _6 s (see Eci l71|l . which is long compared to the correlation time for the vibra- 
tional modes reservoir of 10 _7 s (given by 27r /v max ), and much longer than the 
correlations times associated with the SE or cavity decay modes reservoirs. Some 
relaxation elements obtained are given in Table 4. 



Table 4. Relaxation elements in two qubit gating process 



■T'ia — ; ib— 


r^a — ;jb— ~ (J-'ja— \ ib— ) 


Mo-; jb— 


r C+;C+ 


rfe a -C+: kb-C+ 


4.10 2 s~ i 




4.10 2 s- i 


1.10V 1 


s- 1 










l.lO^s- 1 



The scaling of the decoherence rate with qubit number for tne case of the two 
qubit gating process treated in Tregenna et al gives an overall decoherence 
time which is independent of N. This is because the terms for gated qubits i, j do 
not scale with N, whilst the non-gated qubit contributions (which scale with N) 
are zero for the present case where the 2-1 transition resonantly coupled to cavity 
mode and the 2-0 transition is uncoupled. However, other parameter choices, 
such as having both optical transitions coupled to the cavity mode, could lead 
to large non-gated contributions, due to Tk a -c+; kb-c+ ~ 1.10 +6 s _1 associated 
with LD coupling of the qubits with the cavity and vibrational modes. If one 
photon is present, so that (a^b^b) ~ 1, then even modest size qubit numbers 
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TV ~ 10 3 would lead to non-gated qubit contributions exceeding those from 
other contributions, such as (b^b)Tc+; c+> where Tc+;C+ ~ 1.10 +8 s -1 . 

As for one qubit gating, the contributions from the various terms to (1/ Td) ga ti 
(and hence to the change in fidelity during the two qubit gating period) are sig- 
nificantly different in size. The relaxation element factor may be more important 
than the qubit state factor, and vice versa. We would need to investigate all 
120 terms to determine which contribution is dominent. For example, for the 
terms numbered 23, 24, 34 and 35, the relaxation element is ~ 4.10 2 s~ 1 and the 
qubit factors are ~ 1, giving a product ~ 4.10 s . For term number 120, the 
relaxation element is ~ 1.10 8 s _1 and the qubit factor gives the probability (b^b) 
of finding one photon in the cavity mode. If this probability is greater than 
~ 4.10 -6 (and it could be as high as unity if decoherence free subspaces were 
not utilised during the gating process), the term 120 would be more important 
than the terms 23, 24, 34 and 35. Term 120 is due to cavity decay. If the cavity 
decay term was the most important, the reduction in fidelity would be given by 

AF = -((tfb)-(tf)(b))T cavitv ^- (72) 

which is ~ 6.10 2 (b'b). The probability of finding a photon in the cavity mode 
must be less than 10 -5 if the fidelity loss is to be reasonably small. 



4 Discussion 

The scaling of decoherence effects in circuit model quantum computers have 
been studied for the situation where the number of qubits N becomes large. 
Decoherence effects were specified via the fidelity, with its initial rate of change 
defining the decoherence time scale. Expressions for the decoherence time scale 
were obtained for the intermediate time regime via Markovian theory. The 
general case was treated where the qubit system was in any pure state, the 
reservoirs being in thermal states. The decoherence time scale was expressed 
in terms of Markovian relaxation elements and expectation values of products 
of fluctuation operators for the decohering quantum system. The expression 
given in EqEOI for the decoherence time scale is quite general and may have 
applications for treating decoherence in other macroscopic systems, such as 
Bose condensates or superconductors or in quantum measurement theory. 

A standard model involving N ionic qubits, each a three-state lambda sys- 
tem, was studied, with localised, well-separated qubits undergoing vibrational 
motion in a lattice of trapping potentials. Coherent one and two qubit gating 
processes were controlled by time dependent localised classical EM fields, the 
two qubit gating processes being facilitated by a high Q cavity mode. The qubits 
were coupled to reservoir of spontaneous emission (SE) modes, the cavity mode 
was coupled to a bath of cavity decay modes. For ionic qubits, the numerous 
collective vibrational qubit modes also acted as a reservoir, with Lamb-Dicke 



20 



coupling to the internal qubit system. A key objective of the work was to in- 
vestigate decoherence effects due to the qubit vibrational motion. Parameters 
similar to those in the model treated by Tregenna et al [SS] were chosen, with 
comparable cavity decay and cavity Rabi frequencies, both much larger than 
the spontaneous emission decay rate and the Rabi frequencies of the two qubit 
gating fields. One optical transition was resonant with the cavity mode. For 
two qubit gating, the other transition was also assumed not coupled to the cav- 
ity mode. Our primary aim was to evaluate fundamental rather than technical 
causes of decoherence in standard qubit based quantum computers. 

For the standard model we investigated, cavity decay, spontaneous emission 
and Lamb-Dicke coupling to the vibrational modes were the most important 
fundamental causes of decoherence. Effects due to Rontgen and diamagnetic 
interactions were found to be negligible. Technical causes of decoherence, such 
as fluctuations in the trapping fields, though needed to relate decoherence times 
to current experiments |59| were not studied here. 

Characteristic decoherence time scales were evaluated for specific qubit states 
(Hadamard, GHZ) at finite temperature in the situation with no gating pro- 
cesses occurring. The decoherence time scaled as 1/iV. The decoherence time 
scale for the uncorrelated Hadamard state could be made infinite by choosing 
two optical dipole matrix elements that added to zero. The decoherence time 
scale for the GHZ state was very long, due to the Boltzmann factor - td being 
about 10 19 s for N wlO 4 qubits, even if the free SE decay rate of 10 8 s -1 applied. 

For the case of one qubit gating processes due to weak two photon resonant 
Raman fields with a large one photon detuning, the decoherence time scale was 
evaluated, but at zero temperature. Decoherence was mainly due to Lamb-Dicke 
coupling of the gated qubit with the Raman fields, but the loss of fidelity during 
the gating process was small, being proportional to the square of the Lamb- 
Dicke parameter. Scaling effects were associated with non-gated qubits and were 
small. For both optical transitions coupled to the cavity mode, decoherence was 
associated with Lamb-Dicke coupling of non-gated qubits with the cavity mode, 
no photon being present. However, the effects were negligible even for N wlO 4 
qubits. Scaling effects were absent for only one coupled transition. 

For the case of two qubit gating processes due to one photon resonant gating 
fields, as in the work of Tregenna et al [SSjj the decoherence time scale was eval- 
uated, also at zero temperature. Scaling effects were absent for the parameters 
chosen, so overall the decoherence time is independent of qubit numbers. How- 
ever, other parameter choices, such as having both optical transitions coupled to 
the cavity mode would lead to a significant contribution associated with Lamb- 
Dicke coupling with the cavity mode, one photon being present. In this case, 
modest qubit numbers N wlO 3 qubits would result in non-gated contributions 
that exceed those for the gated qubits. 

In our model the parameters used have been the same as for the theoretical 
model studied by Tregenna et al |5fi| . rather than those where real ions are 
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involved. This was done in order to compare for quantum computer models of 
this type, the effects of including (or not including) the scaling up of qubit num- 
bers and allowing for decoherence due to vibrational motion. A treatment for 
real ions based on three state lambda systems and involving only electric dipole 
transitions is generally too simplified. The presence of other states (such as 
additional magnetic substates, or states associated with other hyperfine levels) 
may need to be taken into account, the actual transitions involved may be of 
electric quadrupole or magnetic dipole character, and magnetic fields may need 
to be present in order that only the 0-2 and 1-2 transitions are resonant with 
the two-qubit gating laser fields. There are several possibilities which involve 
storing the qubit in states 0, 1 and utilising an excited state 2 in the gating pro- 
cesses, so that although these key states form a lambda system, other states or 
non electric dipole transitions may be involved. Consider the case where states 
0, 1 are associated with two hyperfine sublevels of a ground electronic level and 
state 2 is an optical excited state. A simple system of this type involves a 2 Si/ 2 
ground level and a 2 P°/ 2 excited state, but with a non-zero nuclear spin 1 = 1/2. 
With 2 S 1/2 (F = 0, Mp = 0) as state 0, 2 S 1/2 (F = 1,M F = +1) as state 1 and 
with 2 PJ/ 2 (F = l,Mp = +1) as state 2, suitable polarisations for the gating 
laser fields can be chosen to only cause transitions between these states. How- 
ever, spontaneous emission causes transitions into the 2 S!/ 2 (F = l,Mp = 0) 
state, so there is no longer a three state lambda system. A second example is 
where states 0, 1 are associated with two Zeeman substates of a ground elec- 
tronic level and state 2 is an optical excited state. A simple system of this type 
exists in 40 Ca + where states 0, 1 are the ground level 2 Si/ 2 (Mj = —1/2) and 
(Mj = +1/2) states and state 2 is say the metastable 2 D 5 / 2 (Mj = +1/2) state. 
Here the nuclear spin is I = 0, so no hyperfine structure is involved. With suit- 
able polarisations the required one and two qubit gating processes that do not 
involve other states can be performed, the presence of a non-zero magnetic field 
detuning the transition between 2 Si/ 2 (Mj = —1/2) and the additional 2 D 5 / 2 
(Mj = —1/2) state. Spontaneous emission from state 2 only causes transitions 
to states 0, 1, so here a genuine lambda system is involved. However, an elec- 
tric quadrupole transition connects state 2 with and 1 rather than an electric 
dipole transition. If the state 2 was chosen as say the lowest 2 Pg/ 2 (Mj = +1/2) 
state so that electric dipole transitions are involved, then spontaneous emission 
processes to 2 D 5 / 2 and 2 D 3 / 2 states occur and more than three states would be 
involved. A final example involves storing the qubit in states 0, 1 where is a 
ground state and state 1 is a metastable excited state, so that lo\q is an optical 
rather than a Zeeman or hyperfine frequency. Such a case exists in 40 Ca + with 
the states 0, 1 and 2 being magnetic substates of the lowest 2 Si/ 2 , 2 D 5 / 2 and 
2 Pgy 2 energy levels (with respective substates Mj = 1/2,5/2,3/2 for example). 
However, even with suitable laser field polarisations for the one and two qubit 
gating fields, the additional Mj = 1/2, 3/2 substates of the 2 D 5 / 2 level become 
involved due to spontaneous emission from the 2 Py 2 = ^/2) state. Thus, 
the theory would need to be extended to allow for the actual states and radiative 
transitions involved for a particular ion of interest. This choice of ion would be 
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made to minimise the numbers of states needed - suggesting avoidance of cases 
where there are many lower energy (fine structure, hypcrfine stucture) levels, 
together with as small an upper state spontaneous decay rate as possible - sug- 
gesting avoidance of electric dipole downward transitions in favour of electric 
quadrupolc or magnetic dipole processes. The 0-2, 1-2 transitions also need to 
be in the optical frequency range in order to couple these transitions to a high 
Q cavity. Also, cases where there are other levels between 0,1 and 2 are also 
unfavourable, as other such states may be populated via spontaneous emission. 

In conclusion, lambda systems localised in a high Q cavity, which can both 
facilitate two qubit gating processes and reduce decoherence caused by sponta- 
neous emission and cavity decay, are a useful system for research on scalable 
quantum computers. However, for real ions the model needs to be expanded to 
include the presence of all magnetic substates and to treat the case of electric 
quadrupole and magnetic dipole transitions. The case of neutral qubits, where 
the vibrational modes are independent and do not constitute a reservoir, is also 
of significant interest and a treatment via the present Markovian theory would 
be worthwhile. 
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6 Figure caption 

Figure 1. Model of an N qubit quantum computer. Three-state lambda system 
qubits are localised around well-separated positions via trapping potentials, and 
undergo collective centre of mass (CM) vibrational motions. Coherent one and 
two qubit gating processes are controlled by time dependent localised classical 
electromagnetic (EM) fields that address specific qubits. Two qubit gating pro- 
cesses are facilitated by a cavity mode ancilla, which permits state interchange 
between qubits. The lambda system qubits are coupled to a bath of EM field 
spontaneous emission (SE) modes, and the cavity mode is coupled to a bath of 
cavity decay modes. For large N the numerous collective vibrational modes of 
the qubits also act as a reservoir, coupled to the qubits, the cavity mode and 
the SE modes. 
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